AD-A070  200 


UNCLASSIFIED 


WISCONSIN  UN I V— MAD I SON  MATHEMATICS  RESEARCH  CENTER  F/6  20/11 

APEX  SINGULARITIES  FOR  CORNER  CRACKS  UNDER  OPENING*  SLIDING*  AN — ETC(U) 
MAR  79  B NOBLE*  M A HUSSIAN*  S L PU  DAA629-75-C-0024 

MRC-TSR-1940  mi 


BDC  FILE  COPYj 


Mathematics  Research  Center 
University  of  Wisconsin— Madison 

610  Walnut  Street 
Madison,  Wisconsin  53706 


March  1979 


(Received  February  5,  1979) 


D D C 


B 

Approved  for  public  release 
Distribution  unlimited 


Sponsored  by 

U.S.  Army  Research  Office 


P.O.  Box  12211  m 

Research  Triangle  Park 
North  Carolina  27709 


UNIVERSITY  OF  WISCONSIN  - MADISON 
MATHEMATICS  RESEARCH  CENTER 

APEX  SINGULARITIES  FOR  CORNER  CRACKS  UNDER 
OPENING,  SLIDING,  AND  TEARING  MODES 

t t 

B.  Noble,  M.  A.  Hussain  and  S.  L.  Pu 

Technical  Summary  Report  #1940 
March  1979 

ABSTRACT 

Power  singularities  at  the  apex  of  a flat,  wedge-shaped,  angula 
sector  crack  in  a three-dimensional  solid  are  studied.  Using  Boussiu.ng 
stress  functions  of  elasticity  theory,  the  problem  is  reduced  to  an 
eigenvalue  problem  of  a dual  series.  The  stress  singularity  is  found 
to  be  stronger  or  weaker  than  one  half,  depending  upon  whether  the  apex 
angle  is  greater  or  less  than  180°.  This  tends  to  accelerate  or  retard 
the  crack  growth  at  the  apex  until  the  crack  front  straightens  out. 
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APEX  SINGULARITIES  FOR  CORNER  CRACKS  UNDER 
OPENING,  SLIDING,  AND  TEARING  MODES 


B.  Noble,  M.  A.  Hussain'  and  S.  L.  Pu' 


INTRODUCTION 

In  this  paper,  the  stress  singularities  at  the  apex  of  a flat,  ^ 
shaped,  angular  sector  crack  in  a three-dimensional  solid  are  studied  Hu- 
crack  is  subjected  to  an  opening,  sliding,  or  tearing  mode  loading  con  Mon. 
These  modes,  together  with  their  two-dimensional  idealizations  as  commonly 


seen  in  the  literature,  are  shown  in  Figure  1. 
identical  when  the  apex  angle  approaches  180°. 


The  two  sets  of  figure.--  become 


For  the  two-dimensional  case,  it  is  known  that  the  stress  singulnt  itv  at 
the  crack  tip  is  of  the  order  of  one  half  for  all  of  the  modes.  Howev, , , when 
the  two  crack  fronts  meet  at  a sharp  corner,  as  shown  in  Figure  1,  the  stress 
singularity  is  no  longer  of  the  order  of  one  half  but  depends  mainly  on  the 

Included  angle  as  well  as  the  mode  the  crack  is  subjected  to,  namely  mode  1, 

II,  or  111.  A crack  subjected  to  an  arbitrary  load  can  be  decomposed  into 
these  three  basic  modes. 

As  tar  as  singularities  are  concerned,  it  is  of  interest  to  note  tlu 

equivalence  of  crack  and  punch  problems.  In  such  a case,  the  angle  ot  tin 

base  of  the  punch  and  that  of  the  cracked  area  are  explements  ot  each  ollu-i 
and  in  the  loading  process  the  punch  is  held  stationary.  A mode  I crack 
problem  can  be  reduced  to  a problem  in  potential  theory,  but  the  othei  modi 
cannot . 
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Fig.  (1)  - Three  basic  modes  of  two-dimensional,  and  the  correspond- 
ing three-dimensional  flat  wedge-shaped,  crack  surface 
displacements. 


The  potential  problem  of  charge  distribution  at  the  tip  of  a flat  angular 
sector  was  solved  by  Noble  (1959)  [1].  Subsequently,  due  to  a wide  application 
of  potential  theory  in  the  field  of  mathematical  physics,  a number  of  papers  [2- 
7]  appeared  for  the  same  problem.  Even  though  the  methods  and  approaches  are 
different  in  these  papers,  the  final  results  are  the  same  as  [1], 

In  this  paper,  we  use  Boussincsq  potentials  of  the  three-dimensional  theory 
of  elasticity.  Assuming  power  law  singularity  and  using  separation  of  vari- 
ables, the  problem  is  reduced  to  a dual  series  relation  for  the  case  of  mode  I, 
and  coupled  dual  series  for  modes  II  and  III.  Based  on  [1)  a method  is  devel- 
oped to  solve  these  coupled  dual  scries.  The  power,  versatility,  and  gener- 
ality of  the  method  can  be  seen  from  the  simplicity  with  which  new  results  are 
obtained.  It  is  shown  how  results  of  any  desired  degree  of  accuracy  can  be 
obtained  from  simple  algebraic  computations. 
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HU-'  BOUSSINP.SQ  SOl.UTIONS  IN  Sl'ltl  U1CA1.  COORDINATES 


In  tlio  absence  of  body  forces,  the  equation  of  equilibrium  for  a homoge 
neous,  isotropic,  clastic  body,  in  terms  of  displacement  vector,  is 

V*G  ♦ (l-2v)'1VV’u  - 0 (I) 


Here  V is  Poisson's  ratio.  According  to  Bouss inesq,  the  general  elution  of 
nq.  (1)  may  be  written  as  a superposition  of  three  displacement  fields,  iij, 

u2.  a3 

2Gut  * Vij , , 2Gu2  =>  2VX(kO)  , 2Giij  « V(;X)  - 4(1  v)\k  (2) 

where  G is  the  shear  modulus,  and  <j>,  0,  X are  arbitrary  harmonic  functions. 
These  will  be  referred  to  as  basic  solutions  1,  2,  and  A,  respectively. 


We  wish  to  investigate  the  order  of  power  singularities  at  the  apex.  To 
this  end,  it  is  necessary  to  investigate  the  near  field  solution  in  spherical 
coordinates  (r,0,$l,  without  implying  that  the  solid  under  cons iderat  ion  is  of 
a spherical  shape. 

We  choose  the  following  spherical  harmonics 

cos 


ij1  = r^1**  Pra  . (cost)) 
' m \i  ► 1 


0 - r11*1  Pm  , (cosOl 

m p*  1 


Pm (cosO) 
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Pig.  (2)  A flat,  wedge  shaped  crack. 


where  P™  is  the  associated  l.egeudrc  function  of  the  first  kind.  Upon  subst  i 
tut  ion  i|  )t)  into  the  first  basic  solution  of  (21,  we  obtain  a solution  desig 
nated  symool  ical  ly  by  (1  . Similarly  we  obtain  solutions  (2  Om]  and  |A  \m]  . 

The  selection  of  cos  m^  or  sin  tiit}'  in  each  stress  function  depends  on  the 
symmetry  of  the  problem.  Our  objective  is  to  find  the  value  ot  p in  the  open 
range  of  (0,1).  The  final  solution  (Sj  is  the  superposition  of  these  solutions 

[S]  F {A  [l-t|>  1 ♦ B 12-GJ  ♦ D [A  X ]}  (A) 

11  ‘■  m1  in1  m1  m‘  in' 

where,  and  thereafter,  \ denotes  the  summation  with  respect  to  m for  m 0,1, 


On  the  plane  0 3 it/ 2,  the  components  of  displacement  and  the  pertinent  compo- 
nents of  stress  for  [S]  are: 

2Gur  3 l * 2mBm]rWP"+1 

:r,u0  = l l-(ramM)AB  ♦ (3-4v)DmlrUP™ 

**♦  3 H*  "*»  - 2(Wl)BmlryP;M  ^ * 

> (4) 

a0  3 l [-(m^lJA,,  ♦ 2(l-v)Dm(-m*p*l)]ry~1pJ+1  m<> 

T0r  3 I [-p(m*p+l)Am  ; ♦ (l-2v)yDm] r*'~  1 P™  m# 

/ 

T0$  3 I [tm(m*vi*l)Am  ♦ ♦ (l-2v)mDj  rM~  lpJJ  mp 

where  P®  = P™(0)  and  whenever  there  are  two  signs  preceding  a quantity,  the 
sign  on  the  ^top  goes  with  the  trigonometric  function  on  the  top  and  vice  versa. 

1'HRl.H  MOPES  OF  CRACKS  AND  DUAL  SERIES  RELATIONS 

For  a crack  shown  in  Figure  2,  the  leading  edges  of  the  crack  are  <J>  = t ci, 
and  the  crack  is  in  the  x-y  plane  (0  = tt/2).  Let  D“  and  D+  be  the  cracked  and 
uncracked  region  of  the  plane  0 3 ir/2.  Within  the  cracked  region,  the  displace- 
ment is  discontinuous.  If  the  discontinuity  is  in  the  z-direction  (u0  - u0  3 
finite),  the  crack  is  under  mode  I;  if  the  discontinuity  is  in  the  x-uirection 
(u*  - u~  3 finite),  the  crack  is  defined  to  be  under  mode  II;  and  if  Uy  - Uy  = 
finite,  the  crack  is  defined  to  be  under  mode  III.  Boundary  conditions  for 
various  modes  are  tabulated  below. 


15) 

(6) 

(7) 


Boundary  Conditions  on  0 = v/2 


Non  Mixed  Conditions 
(in  D"  ♦ D*) 

Mixed  Conditions 

in  D" 

in  D* 

Mode  I 

Mode  II 

Mode  III 

T0r  3 T0<|»  3 0 

O0  = 0 

o0  3 0 

oQ  3 0 

T0r  3 Te<J>  = 0 

T0r  3 Te*  = 0 

u0  3 ° 
ur  = u^  = 0 

ur  3 u $ 3 0 

For  mode  I,  ug  is  even  in  $.  This  leads  to  the  use  of  the  trigonometric 
functions  at  the  top  of  (4).  The  boundary  conditions  of  (6)  and  (7)  are 
identical,  but  the  symmetric  properties  are  different  for  mode  II  and  mode 
III.  In  the  former  case,  ur  is  even  and  u^  is  odd  in  $ while  in  the  latter, 
the  reverse  is  true.  Hence  the  proper  set  of  quantities  should  be  selected 
in  (4),  for  each  case. 
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A.  Movie  l 


Using  l'q.  (1)  .tml  the  non  mixed  conditions  of  (S),  we  li.no 
Rm  * 0 and  Am  =•  (imp*  1 ) ' 1 ( 1 -2v)Dm  ( 

The  mixed  bound  ry  conditions  of  (5)  and  using  (8)  and  (4),  yield 
J bm  cosm$  * 0 0<$<a 

l 

l q„b|n  cosing  =*  0 ix<4»<it 

where 

bm  = (-B*y*l)DmPjtl  , qm  - (-m»H*l)  P"/Pjj+1  O 

R.  Mode  II 

For  the  homogeneous  condition  of  Og  in  (6)  vising  (4)  the  »■•«•<*'  • • 

A and  D must  satisfy 

A|n  =*  (ra+M+1)  l2(l-v)Dm 

This  relation  and  the  mixed  conditions  of  (6)  yield  the  following 
se ties : 


) Fmcosm$  3 0 

0<4><a 

l <RmFVSmFm>C0:im*  5 0 

«<^<n 

£ Fmsiim|>  a 0 

0<$<a 

?I(UmW‘»)sinm*  “ 0 

a<4><ir 

where  denotes  the  summation  with  respect  to  m for  m 1,2, 

Fra  * (m*u*l)[liAra*2m(l-v)RjP™  , F)n  ■ (m+y+l)  [m\>*2)i(l-v)Rnll,J‘l  l 

- lm*-n(n>n(l-v)lVm  , Sm  - m(l-V-vp)Vm  1 

Tm  » m(Hyv)Vm  , Um  =>  [(l-v)ra*  - li(u+i)]Vm  0 

in  (14)  and  IIS)  \'m  stands  for  Pj+1/Pj/l  C»*  V)]  • 

C.  Mode  III 

Similar  to  the  preceding  case,  we  have 
J\F.rasinm*  * 0 0<4»<a 

illR»FVSwF»)slnm*  a 0 rt<*iir 


(17) 


J Fmcosm$  = 0 0<$<ot 

l (UmFra+TraEm^cosm^  3 0 a<*** 

SOLUTIONS  OF  SINK  AND  COSINE  SERIES 

For  t he  solution  of  (9),  (11),  and  (12),  (16)  and  (17),  we  need  the  exact 
solutions  of  certain  dual  sine  and  cosine  sc> ies. 

A.  Sine  Series 

Consider  a dual  series  of  the  form 


XQ<  X<1T 


0<_x<xo 


J^sinnx  = 0 

[j  !{  a„sinnx  = f(x) 

A solution  is  given  by  [8] 

R(t)  = - I sec  t d fc  sinuG(u)du 
11  - Jt  { /cost  - cosu 

UlthG(u)  = cot  H.  d fU  sinXx/llfJxidx 

2 du  ' n — r 

o /(cosx  - cosu) 

"’here  xQ 

g(x)  = [an  sinnx,  for  0<x<xQ,  and  am  = - / g(t)  sinmtdt. 

1 x 

For  f(x)  = Cj  sinx,  we  obtain,  using  the  abbreviation  g = cos  xQ, 


-1/2 


C1  x 

g(x)  = — sin  x-  (2cosx  ♦ 1 - g)  (cosx  - 8) 

/2  2 


*•»  - 2 ^m  - e**,.!  ‘ Pm+1  ♦ Pm.  2> ’ Pm  « Pm^ 

For  f(x)  = C2 sin2x,  the  corresponding  results  are 

g(x)  - - /I  C2sin  ^ (cosx-8)  (P+P2 -2cosx+2gcosx-4cos2x) 
an  3 - \ l-2Pn-3+2Ppn-2+(-l^2)Pn.^(l-82)Fn-28PnM>2Pn+2] 
Similarly,  for  f(x)  = C_jsin3x,  we  have 
3C,sin(x/2) 

g(x)  = [16cosVScos2x(l-g)-2cosx(3*2p+p2)-ps-82  + 3p-l] 

2/2  (cosx  -8) 

an  3 M2(Pn.4-Pu3)-2P(Pn.3-Pn+2)-(P2-l){Pn.2>P(Pn.1  + Pn)-Pn+1}l/4 


(IS) 


(19) 


(20) 


(21) 

(22) 

(23) 

(24) 

(25) 

(26) 


H.  fosine  Seri  os 

The  cosine  series  under  consideration  have  the  foim 


V a,,cosnx  0 


x0<x<it 


(27) 


1 

poao  * £ n ancosnx  3 I Pnancosnx» 


CNx<Xr 


A method  of  solution  is  developed  by  Noble  [1]  for  finding  the  lowest  eigen- 
value for  mixed  problems  in  potential  theory.  The  method  uses  a stretching 
transformation  due  to  Schwinger  [9].  bet 


T afp  osnx  - h(x)  for  0<x<xQ 


t hen 


1 r X o Z f V 

a0  J M*)dx,  an  = ~ J h (x)cosnxdx 


2 f °, 

o 


(28) 

(29) 


P<x<x0. 


P0  ' °J  h(u)du+£  ^ j h(u)cosmu  cosmxdu=£pm  “-/  h(u)cosmudu  cosmx,  (30) 


Upon  substitution  from  (29)  into  the  second  equation  of  (27),  we  have,  for 

N 

h(u)du+)  *■-  J h(u)cosmu  cosmxdu 

o 1 m ir  o 

Using  the  stretching  transformation 
cosu  - t ♦ s cos£ 

where  t and  s,  determined  from  u = 0,  xQ  corresponding  to  F,  = 0,n  are 
t = cos2 (xG/ 2) , s = sin2(x0/2) 

A similar  relation  is  assumed  between  x and  f;.  Using  [10],  we  obtain 

),  ^eosmu  cosmx  = - ^£n{ 2 1 cosu-cosx | }- - ?ns/2-*T , *vosm£  cosmt 
l m 2 1 1 m 

bquat ion  (30)  becomes 

n oo  it 

(p0  Cns)  it/  h('0  ;!■'  M.+1  ' “•/  h(u)  cosm f,df,  cosme  = 

o ^ 1 o * 

Y°m  \l  h(l,i  cosnuulr> 

1 o ^ 


(31) 


(32) 


(33) 


cosmx  , 


0<x<tt 


(341 


On  the  right  hand  side  of  (34),  cosmu  and  cosmx  can  be  transformed  to  cosmf, 
and  cosing  in  virtue  of  (31).  bet 


M«0  l cn  cosnC 


0<C<1T 
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t hen 


l I h(u)  ^ 3 c0  • l'  / h(u)  cosmf.dF.  - c„ 


Upon  substitution  from  (36)  into  (34),  we  have  a trigonometric  equation  in 
cosmC.  fquat  ing  the  coefficients  of  like  terms,  a system  of  algebraic  equa- 
tions in  c's  is  obtained.  For  non-trivial  solutions,  the  determinant  of  the 
set  of  equations  must  vanish  which  is  used  to  find  the  eigenvalues. 

Specifically,  if  pn  = 0 for  n > 1 in  the  last  equation  of  (27),  then  cra  0 
for  m > 1 and  for  cQ  ? 0,  we  must  have 

p0  - tns  = 0 (37) 

For  the  case  pj  / 0 and  Pra  3 0,  m 2 then  cm  3 0 for  m > 2 and 

f-p0  ens  ♦ :t2P!  st p j 1 pc0  1 


Similarly,  when  N =»  2 on  the  right  hand  side  of  (27),  c = 0 for  in  > 3 and 
C0.  C|,  Cj  satisfy,  using  ^ for  s2+2t2-l, 


-p0*tns>2t2p1»2ei2p2  stp1  + 4sti1p2  e iP2 

2stpl*2(4st)e1p2  -l>s2p1+(1st)2p2  s2(4st)p2 

2s2eip2  s2(4st)p2  - \ ♦ s\s 


C ! = 0 


1 w 

7 + s p2 


SOLUTIONS  OF  FLAT,  WF.PGE-SHAPED  CRACKS 


A.  Mode  I 


The  second  equation  of  the  dual  cosine  series  (0)  can  be  written  as 

%bo  + l\  i hm  cosm<^  = cosm*  (4t 

Using  the  last  equation,  p.  63  of  [11],  we  have 

>i,  ■ j rc^tj/ro  . («i 

For  a large  m,  qm  -*■  1/m,  [1],  therefore  the  infinite  series  on  the  right  hand 
side  of  equation  (40)  could  be  truncated. 

In  order  to  use  the  solution  of  cosine  series  presented  in  the  preceding 
section,  we  make  the  following  change  of  variables 

$ 3 u-w,  xQ  a ir-ct,  cosm<{>  = (-l)racosmto,  (-l)"b|n  = a^  (42 
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1'ho  dual  cosine  series  (.1))  take  the  form 


ao  * ham  cosmj  = 0 


X0<U><TT 


+ U i am  3 - llm)3m  cosnVi)  0<u)<xo 


As  a first  approximation,  we  drop  out  t be  infip’te  series  on  the  r Vt  id 
of  (43).  Using  (37),  we  have  a transcendental  equation 


COS(J)  = exp(q0/2),  q0  = - \ [r(^)/r(l  ♦ ^)]2 


from  which  y can  he  determined  for  a given  half  angle  a of  the  weige.  or 

the  second  approximation,  only  one  term  is  taken  in  the  series  on  the 
hand  side  of  (43).  Using  (38)  the  value  of  y is  found  as  follows: 


q0-£ns-2t2 (1-qj)  -st(l-qj) 


-2st ( 1 - q j ) 


i-s2  ci-qx) ! 


= 0 (45) 


where  s,t  are  given  by  (32),  (42)  and  qg, 
qj  from  (41). 

Similarly,  the  third  approximation  can 
be  carried  out.  The  results  of  each  approx- 
imation are  shown  in  Figure  3,  showing  a 
remarkably  fast  convergence. 


\ A i 1 l 

Ob 


Fig.  (3)  - y vs.  a/n  f, 
crack . 


rude  I 


B.  Mode  II 


Using  the  property  F™  . (0)/Pm(0)  -*■  (-1)  for  a large  m, 
asymptotically:  * ^ 

Rm  (-!/»)  . Sm  > (-1/m2)  , Tm  > (-1/m2)  , Um  <•  -(l-v)/r 


lhe  second  equations  of  the  coupled  dual  scries  (11)  and  (12)  can  : 
00  00  1 U T 

£ i Fmsinn,4>  = Z((n-  ♦ imv0Fm  + fly  Em]sinm*  ft<4<7T 


RoKo  * | m Emcosm<|)  = £[(Jr  ♦ R,n^)Km+smFmlcosni<<) 

After  the  change  of  variables  similar  to  (42)  with  sinm<}>  = (-1)' 
( l)mF.m  and  Fm*  = (-l)m+*Fm.  The  dual  series  take  the  form 


_q_ 


JjFn/sinmco  = 0,  x0<w<tt 

h “ Fm*sinmul  = Iiclmsin,n“*  0<UJ<xo 
jEm*cosmto  = 0,  x0<lo<it 

1 r»  ^ 

-RoRo*^l  - Em*C0Sma)  = ^lc2mRm*cos,ntJ*  0<0)<xo 


(48) 


(49) 


1.0 


10.5 


where 


•1 


' lin 


= (1-v)  (U_Fm*-TmEm*) +F_*/m  , 
^ J '•mm  mm-*  m 


v - o. 


-1 


c,  = m +R  -S  F */E  * 
2m  m m m m 


(50) 


a/ if-  - 


For  a large  m,  both  and  c?m  1/m 


'2m 


Fig-  (4)  - y vs-  cl/tt  for 

a mode  II  or  a 
mode  III  crack. 


As  a first  approximation,  we  take  only  one  term,  m=l,  in  the  infinite 
series  on  the  right  hand  sides  of  (48)  and  (49).  From  (22)  we  obtain 

Fx*  = (1-B)(3+B)[F1*-(1-v)'1(U1F1*+T1E1*)]/4 

where  E^*  can  be  determined  from  (29),  (31)  and  (36). 

')  x0  p TT  , 

E^*  = — J h(u)cosudu  = — / h(u)  (t+scos?)d  = 2tCg+sc^ 

77  o o c s 

For  the  cosine  series  (49),  we  obtain  two  equations  from  (45) 

(-R  -£ns-2c_1t2)cri-stc-,c1  = 0 
o 21  ' 0 21  1 

-2stc-?1  cn+  (1-s2c2]  )Ci  = 0 


(51) 


(52) 


These  two  equations  and  (51)  constitute  a system  of  linear  equations  in  eg,  cj, 
-1 

RQ+£ns+2t2  (1  + R-^)  st(l  + R1) 


Fj*.  For  a non-trivial  solution,  we  must  have 


-tSi 


-l+sz (1+Ri ) 


-SSi 


= 0 (53) 


2st(l+Rx) 

-tT1/2  -sTj/4  -(l-v)/(l-3)/(3+e)+(l-v-U1)/4 

This  is  the  equation  for  the  determination  of  y. 

We  can  proceed  in  a similar  manner  to  higher  order  approximations.  In 
Figure  4,  we  have  shown  results  for  v = 0 for  the  first  three  approximations. 
The  difference  between  the  second  and  third  approximations  is  very  small. 
Hence  we  stopped  at  the  third.  Since  the  results  did  not  differ  from  the 
second  to  the  third  approximation,  we  only  plotted  the  third  approximation 
jn  Figure  5 for  different  Poisson's  ratios. 


Fig.  (5)  - The  third  approximation  of  p vs.  a/v  for  a inode  II  or  a mode 
III  crack. 

C . Mode  1 1 1 

In  comparing  equations  (16),  (17)  with  (11)  and  (12),  we  see  that  the 
coupled  dual  series  for  mode  III  are  obtained  from  those  for  mode  11  by  inter 
changing  Em  and  Fm,  Rm  and  U,n,  Sm  and  Tm.  The  same  approximate  method  for 
mode  II  can  be  applied  to  mode  III.  For  v = 0,  mode  II  and  mode  III  have  the 
same  results  which  are  plotted  in  Figure  4.  For  v = 0.25  and  v * 0.5,  only 
results  from  the  third  approximation  are  plotted  in  Figure  5. 

I'ONi : I US  IONS 

The  equivalence  between  the  potential  problem  and  the  crack  problem  under 
mode  I or  the  frictionless  punch  problem  has  been  known  for  a long  time.  The 
wide  applicability  of  potential  problems  has  been  a prime  motivation  in  obtain 
ing  the  singularities  at  the  apex  by  various  methods.  In  this  paper,  however, 
we  have  presented  a method  of  solution  for  mode  II  and  III  conditions.  The 
method  gives  results  to  any  desired  accuracy  and,  to  the  best  of  our  knowl 
edge,  is  not  yet  available  in  the  open  literature. 

For  modes  other  than  opening  mode,  the  results  show  that  the  stress 
singularities  are  dominated  by  the  angle  of  the  apex  as  well  as  the  elastic 
constant  of  t lie  material,  flic  results  further  indicate  that  when  the  apex 
angle  is  greater  than  180°,  the  stress  singularity  is  stronger  than  one  half 
enhancing  the  tendency  of  crack  front  to  straighten  out.  Similarly,  when  the 
apex  angle  is  less  than  180°,  the  stress  singularity  is  less  severe  than  one 
half  and,  again,  this  will  tend  to  retard  the  growth  at  the  apex  until  the 
track  front  straightens  out. 


1 


i 


It  should  be  noted  that  the  definition  of  modes  II  and  111  is  quite 
arbitrary  for  the  configuration  under  consideration,  since  the  crack  fronts 
away  from  the  apex  are  under  mixed  mode  conditions  of  the  conventional  type. 
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